Neutrino spin relaxation in medium with stochastic characteristics 
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The helicity evolution of a neutrino propagating in randomly moving and polarized matter is 
studied. The type of the neutrino interaction with background fermions is arbitrary. We derive the 
equation for the description of the averaged neutrino helicity evolution. In the particular case of a 
r-neutrino interacting with ultrarelativistic electron-positron plasma we obtain the expression for 
the neutrino helicity relaxation rate in the explicit form. We study the neutrino spin relaxation 
in the relativistic primordial plasma. Supposing that the conversion of left-handed neutrinos into 
right-handed ones is suppressed at the early stages of the Universe evolution we get the upper limit 
on the r-neutrino mass. 
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It is well known that if a particle, possessing the mag- 
netic moment, interacts with the external electromag- 
netic field, particle's spin will rotate around a certain 
direction determined by this external field. The covari- 
ant description of the particle's spin precession in the 
external electromagnetic field was presented in Ref. 0. 
Recently it was established (see Refs. 0,0) that weak 
interactions can also cause the particle's spin preces- 
sion. A neutral 1/2 spin particle (e.g., a neutrino) in- 
teracting with matter via weak currents was studied in 
Refs. EL 0- R was found that the particle's spin 

precession rate depends on the velocities and polariza- 
tions of the background matter fermions. On the basis of 
the elaborated approach one can study various interest- 
ing phenomena such as, for instance, neutrino spin light 
and self-polarization effect in matter (see Refs. [a,LZ|)- 

In case of a neutrino interacting with medium the effect 
of velocities and polarizations of the background fermions 
on the neutrino's spin rotation is suppressed by the factor 
I/7, where 7 = E v /m v . Thus for relativistic neutrinos 
the transitions between left and right polarized particles 
due to the matter effects vanish. In order to achieve 
the essential neutrino spin rotation rate one should have 
background fermions with very high velocities or polar- 
izations. Such ultrarelativistic matter can exist mainly in 
various astrophysical environments or in the relativistic 
plasma of the early Universe. 

Neutrinos in the early Universe were discussed for the 
first time in Ref. 0. Since then a great deal of papers 
devoted to this topic have been published. The compre- 
hensive review of neutrinos in cosmology was presented 
in Ref. . The neutrino mass constraints can also be ob- 
tained from the studying of the early Universe evolution. 
For example, one extracts information about neutrino 
mass from the Cosmic Microwave Background radiation, 
power spectrum in large scale structure surveys as well 



as Lyman a forest studies. The review of these issues 
was presented in Ref. [Io| . 

In this paper we study neutrino spin dynamics in 
medium with stochastic characteristics such as num- 
ber density, velocities and polarizations of background 
fermions. The analysis is based on the generalized 
relativistic invariant Bargmann-Michel-Telegdi equation 
which describes neutrino spin evolution in presence of 
general interactions. That is why the type of the neu- 
trino interactions with medium remains arbitrary. This 
equation enables one to study the neutrino spin preces- 
sion at arbitrary moving and polarized matter even at the 
absence of the external electromagnetic fields. This case 
can be of great importance if a neutrino is treated, for 
example, within the context of minimally extended stan- 
dard model. The magnetic moment of such a neutrino 
has quite small value (see, e.g., Refs. PE3), 
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The external electromagnetic field cannot effectively in- 
fluence on the spin dynamics of a neutrino with such 
electromagnetic properties. We derive the most general 
equation for the description of the averaged neutrino he- 
licity evolution. It is found that the neutrino helicity re- 
laxation rate depends on the correlation functions of the 
hydrodynamical currents of background fermions. These 
correlators are calculated for the ulrarelativistic electron- 
positron plasma. The results of the paper are used for the 
evaluation of cosmological upper bound on the r-neutrino 
mass. 

We start with the generalized Bargmann-Michel- 
Telegdi equation for the neutral 1/2 spin fermion (e.g., 
neutrino) interacting with arbitrary moving and polar- 
ized matter (see Ref. 0), 
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where G^ v = (— P, M) is the antisymmetric tensor which 
accounts for the matter effects, G^ u = (\/2)£ llva pG al} — 
(M,P) is the tensor dual to G^, p, and e are the mag- 
netic and electric dipole moments of a particle, is the 
four velocity and is the four dimensional spin vec- 
tor. The derivative is taken with respect to the particle's 
proper time r. 

According to the results of Ref. |J 
be represented in the following way 



the tensor G^ v can 



G tiu 



x/39 



where the four vectors — (<?°,g) and — (/°,f) 
depend on the four vectors of hydrodynamical currents 
and polarizations of medium fermions, 



E 



The constants, p* 1 ' 2 ^ and $~ , are fixed by the model of 
the neutrino interaction with background fermions, the 
sum is taken over all fermionic species. It is also possible 
to represent the dependence of the vectors M and P on 
the components of the vectors <? M and (see Ref. Q), 

M = 7 {(.9°/3-g)-[/3xf]}, 
P = -7{(/°/3-f) + [/3xg]}, 

where (3 is the neutrino velocity. The dependence of 
these vectors on the velocities and polarizations of the 
background fermions is also presented in Ref. |4| in the 
explicit form. 

It is convenient to rescale the vectors M and P to ex- 
clude magnetic and electric dipole moments from Eq. J5J , 



M 



M 

/i 



P 



Then we should rewrite Eq. using three dimensional 
neutrino spin vector, which is related to the four di- 
mensional neutrino spin vector by the following formula 



p(Cp) 



\ m v ' m v (m v + EJ) J 1 



(3) 



where p is the neutrino momentum. With help of Eqs. (0) 
and J3} we obtain the equation for the vector 



§ = -[CxM ]-^[CxP ] 

at 7 7 



(4) 



Note that the dynamics of the neutrino spin is deter- 
mined by the vectors M and P in the neutrino rest frame, 



M = 7/3 9 



1 



1 + 7" 



n(g/3) -g 



Po = -7/3 / 



1 + 7" 



It is convenient to introduce the new vector Jl = Mo — Po 
and rewrite Eq. (@J in the form, 



dt 



7 



-[Cxfi]. 



(5) 



The neutrino spin precession is now determined by the 
vector f2 which can be expressed in terms of the com- 
ponents of the new auxiliary four vector — (A , A), 



fi = 7/3 [A - 



1 + 7" 



-(A/3) -A 



(0) 



The four vector A^ also depends on the four vectors of 
hydrodynamical currents and polarizations of medium 
fermions, 



where the coefficients pf and £/ are related to the con- 
stants Py 1 ' 2 ^ and £f\ 

Pf - Pf + Pf , ?/ - ?/ +?/ • 

Now let us turn to the description of the neutrino helic- 
ity evolution. Basing upon Eqs. © and © one obtains 
the equation for particle's hclicity, h = (f3Q/f3, which 
has the form 

^ = ^n[Cxfl]) = ~(A[nxC]), (7) 

where n = (3/(3 is the unit vector along the neutrino 
velocity. 

We can solve Eq. @ for the description of the neutrino 
helicity evolution by means of iterations. The similar 
method for the description of the neutrino helicity evo- 
lution in stochastic electromagnetic fields was proposed 
in Ref. [R|. One substitutes the formal integral solution 
of Eq. ©, 

C(t) = - [ dt'[C(t') x n(t% 

7 Jo 

in Eq. J7J. After some uncomplicated but rather cumber- 
some calculations we obtain the equation for the neutrino 
helicity, 



dt'{(A ± (r,t)A ± (r',t'))h(t') 



dh _ 

dt \7/ Jo 

+ 7 (Cx(f)A ± (r,t)) [/3^ (r',O-A||(r',t')] }• (8) 



The longitudinal, A\\ = (nA), and transversal, A^ = 
A — n(nA), components of the vector A are taken with 
respect to the neutrino velocity. The vectors r and r' 
are the neutrino positions at the times t and t' . If we 
suppose that a neutrino propagates along a line, then 
r' - r = 0(t' - t). 
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Now let us discuss the average quantities in Eq. JSJ 
that will be denoted as (...). We have already men- 
tioned that the neutrino helicity evolution is suppressed 
by the factor I/7 which is small for relativistic neutrinos. 
This feature can be also derived from Eq. JSJ. Thus, if 
we consider Eq. JHJ to the lowest order in 1/7, the vari- 
ation of the averaged neutrino helicity turns out to be 
negligibly small during the correlation time of the vector 
components, because the characteristic time scale of 
the correlators (A M (r, t)A v (r' , t')) for the high tempera- 
ture plasma is very short. Let us evaluate this time scale 
of the correlators when a r-neutrino interacts with high 
temperature electron-positron plasma (this problem will 
be discussed in details below). The time scale is approx- 
imately equal to the mean free path of plasma particles 
(see, e.g., Ref. Q)i i- e - electrons and positrons, 



t 



T 2 



where T is the plasma temperature, n in the number 
density in plasma, a cm = e 2 is the fine structure constant 
and 



A = In 



T 



T 



4ira cm n 



is the Coulomb logarithm. 

Supposing that for the relativistic plasma one has 
n pa 0.36T 3 (here we take into account that both elec- 
trons and positrons conribute to the number density) 
and setting T ps 100 MeV, we get for the mean free 
path the following estimate, t « 5.2 x 10~ 20 s. On the 
other hand one obtains for the typical time of the av- 
eraged neutrino helicity relaxation (see Eq. l(T7|) below), 
1/r pa 8.6 x 10 -4 s. One can see from this estimates that 
t <C 1/r. That is why we can consider the integrand in 
Eq. (JSJ) as the product of two functions: the correlators 
of the vector components and the averaged neutrino 
helicity. 

The vector A M is the linear combination of the hydro- 
dynamical currents and polarizations vectors. Different 
components of these vectors are sure to be uncorrelated. 
Thus we get that the correlators of different components 
of the vector A 11 vanish. Finally one obtains more simple 
equation for the averaged neutrino helicity 



d(h) 
dt 



dt'{A ± (r,t)A ± (r',t')){h(t')). (9) 



Eq. © is the most general one which describes the av- 
eraged neutrino helicity evolution in presence of mat- 
ter with random velocities and polarizations. Note that 
the type of the neutrino interaction with background 
fcrmions is not fixed yet. 

Eq. @ can be transformed to the more convenient 
form which will be used in the further analysis. As we 
have already noticed the correlations of the fields A are 
not equal to zero for a rather short period of time. Thus 
supposing again that the neutrino helicity is practically 



constant during the matter fields correlation period as 
well as taking into account time-translation invariance 
of the matter fields correlators we receive more simple 
differential equation for the mean neutrino helicity 



d(h) 
dt 



= -r(h), 



where neutrino helicity relaxation rate has the form 



di(A ± (r,t)A ± (Q,0))| r=/9t . (10) 
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We can take into account only (j/(r, f)j/(0, 0)) correla- 
tion functions in Eq. (|10fl if we consider neutrino inter- 
acting with the relativistic plasma. Indeed for relativis- 
tic fermions the polarization vector is proportional to the 
current, A^ w /i/j/, where hf = (/3/C/)//?/ is the fermion 
helicity, j3f is the velocity of the reference frame in which 
the mean momentum of /-type background fermions is 
zero and £/ is the mean value of the polarization vector of 
the /-type background fermions in the above mentioned 
reference frame. The fermion helicities and currents in a 
stochastic plasma take the form 

hf = (h f )+Sh f) jf = (j f ) +6j f , 

where (hf) — and (jf) = are the mean values of 
the corresponding quantities. The fluctuations of fermion 
helicities, Shf, and currents, 5jf, are taken to be small. 
Therefore the correlation functions like (A/(r, <)j/(0, 0)) 
and (A/(r, t)A/(0, 0)) are next to the leading order. 

The hydrodynamical currents correlation functions can 
be expressed in terms of the electromagnetic currents cor- 
relators for plasma consisted of charged particles. If the 
charges of particles in plasma are equal to ±e, then 

(j/(r,t)j/(0,0)> = itirM)j/ m ( > )>> 

where j" n (r, t) is the electromagnetic current. Thus the 
neutrino helicity relaxation rate is expressed in the fol- 
lowing way 



x / d<a^(r,t)j}1(0,0))| p . si . (11) 



The correlation functions in Eq. 111(1 are usually ex- 
pressed in terms of the Fourier transforms, 

(ir(M)jr(°>°)> 

= 7^ J dujd 3 k(jrj- m u^ Mkr - (12) 

The electromagnetic currents correlators can be obtained 
with help of the fluctuation-dissipation theorem (see 
Ref. ^H)- The explicit form of the correlation functions 
is presented in Ref. [T^ . 
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(jj jj )wk 



1 



27rexp(w/T) - 1 



k 2 \e,\ 2 



~fc2~ 



{i-ik/uff 



\e t -(k/co) 2 \ 2 



(13) 



The correlators are expressed in terms of the permittivity tensor components, ei and e t . For the case of the ultrarela- 
tivistic plasma, T>mj, where m/ is the background fermion mass, the longitudinal and transversal components of 
the permittivity tensor take the form (see, e.g., Ref. |l7jl 



e*(k,w) 



1 







~\k-w\~ 










2A: 


k + LU 



2a 2 k 2 



4fc 



1 _ 2ZJ (1 " (Lj/fc)2) ln 



7TW . . 

■ + sign(/c — uj)\ 



irk 



(14) 



- %— (1 - Mfc) 2 ) [1 + sign(/c - W )] [> , (15) 



where a is the Debye length which is related to the par- 
ticles densities in plasma by the following formula 



and 



sign(x) 



47re 2 



f 




is the step function. 

As it was pointed out in Ref. [l3T |. we can take into 
account only longitudinal component of the electromag- 
netic currents correlators since the imaginary part of the 
transversal component of the permittivity tensor [which 
contributes to the transversal component of the electro- 
magnetic currents correlators, see Eq. I|13[) ] is suppressed 
by the additional small factor I/7 2 . Using Eqs. i|ll[l - H15|) 
one obtains for the neutrino helicity relaxation rate, 



iViV 2 t 3 . f +1 

7 J e 2 ^ Pf 16n 2 a 2 J _^ 



f°° LU 1 

d(cos i?) sin 2 cos / dfc - — — -r— 

Jo exp(w/T) - 1 \ei\ 2 



2' 

1J • ' 



uj—k(3 cos i5 



(16) 



where N is the number of particles species in plasma and 
■d is the angle between vectors k and (3. In Eq. (|16fl we 
drop small terms like 1 /j 4 as well as we keep the leading 
terms in the fine structure constant. We also suppose the 
frequency to and the wave vector k are connected by the 
relation, u> = (k/3) = fc/3cosi9, that immediately follows 
from Eq. (fTTfl . 



The neutrino helicity relaxation due to the neutrino 
magnetic moment interaction with the stochastic electro- 
magnetic field was studied in Ref. The helicity re- 
laxation rate obtained in that paper is proportional to /i 2 , 
where [i is the neutrino magnetic moment. One can also 
consider the correlations between stochastic electromag- 
netic fields and polarizations of the background fermions. 
Indeed the background fermions are polarized under the 
influence of the external electromagnetic field. In this 
case the induced polarization of the /-type fermions has 



the form, 



/•(ii 



-d) 




Xf 



(m) 1 
9 f 



-B 



Xf 



where Kf and \f are dielectric and magnetic susceptibil- 
ities, gf' and are the gyroelectric and gyromagnetic 
ratios, rif is the number density of the /-type background 
fermions. In the latter case the corresponding contribu- 
tion to the neutrino helicity relaxation rate will be pro- 
portional to the factor /1 x 7 -1 . 

Supposing that neutrino magnetic moment is described 
within a realistic model, for instance the minimally ex- 
tended standard model (see Refs. 0, 0|), we obtain, 
e.g., for a r- neutrino (this neutrino type will be con- 
sidered below) magnetic moment fi Vr ps 5.8 x 10~ 12 /is 
[see also Eq. Q]. The r-neutrino mass is taken to be 
18.2 MeV [l8| . If a particle has such a small magnetic mo- 
ment, the contribution of the corresponding terms to the 
neutrino helicity relaxation rate will be negligibly small 
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compared to the contribution examined in our paper. 

To evaluate the mean neutrino energy in Eq. Ijl6(l we 
note that it can be expressed with help of the neutrino 
energy and number densities {E v ) — QE/n v . For the 
relativistic neutrino gas which is in equilibrium at the 
temperature T we have 



QE = 2 



d 3 p 



and 



(2tt) 3 exp(£„/T) + 1 
d 3 p 



1 



(2tt) 3 exp(E„/T) + 1 



0.58T 4 



0.18T 3 



(17) 



(18) 



Finally we get for the mean neutrino energy (E u ) ~ 
3.22T. 

Let us discuss a r-neutrino propagating in ultrarel- 
ativistic plasma composed of electrons and positrons. 
This neutrino flavor eigenstate interacts with electron- 
positron plasma by means of the weak neutral currents. 
In this case the coefficients pf are expressed in the fol- 
lowing way 



-Pe+ 



G p 
2V2 



(1 -4sin 2 6» w ), (19) 



where Gf is the Fermi constant and 9w is the Weinberg 
angle. Using Eqs. lfll)jl - lfT§|) we obtain for the neutrino 
helicity relaxation rate, 



T ps 1.16 x 10' 



V 10 Me 1 



T 



MeV J V 100 MeV 



(20) 



The main results of the paper can be applied, e.g., 
for the description of the neutrino spin relaxation in the 
primordial plasma. The number of right-handed neu- 
trinos is very small at present. Indeed if one observed 
a right-handed neutrino, it would be a direct indica- 
tion on the neutrino having a non-zero mass since the 
wave function of a massive particle can correspond to 
both left- and right-handed helicity states contrary to 
the case of a massless particle. However, up to now only 
the upper bound on the neutrino mass is experimentally 
established. Thus we should assert that if there were 
right-handed neutrinos, their number would be not very 
great. The neutrinos with this polarization cannot be 
produced in great amount at the early stages of the Uni- 
verse evolution. Thus we should suppose that the pro- 
cess of neutrino helicity relaxation (i.e. the production of 
right-handed neutrinos since all neutrinos are taken to be 
left-handed initially) is out of the thermodynamical equi- 
librium, i.e. r < H, where H is the Hubble constant. At 
radiation dominated era of the early Universe the Hub- 
ble constant can be approximately evaluated (see, for in- 
stance, Ref. 



H w (G N ) 1/2 T 2 re 1.24 x 10 3 



T 



100 MeV 



(21) 



where Gat is the Newton's constant. 

When temperature of the early Universe is equal or less 
than 100 MeV, the muons and antimuons densities are 
negligible to effectively contribute to the neutrino helic- 
ity relaxation. The densities of neutrons and protons are 
also small. Thus a r-neutrino can interact only with elec- 
trons and positrons at this temperature. From Eqs. I|20|) 
and (ED for T = 100 MeV we get < 10.2 MeV. The 
obtained upper limit on the r-neutrino mass is compat- 
ible with the up to date neutrino mass bounds derived 
in various techniques (see Ref. 0). It should be noted 
that the calculations made in this paper are also consis- 
tent with the obtained upper limit. Indeed in deriving 
of Eq. JTHJ we supposed that I/7 4 = {m u /E v ) A < 1. 
Eqs. (O and JTSJ are valid while m„ < T. All these 
assumptions are in agreement with the temperature of 
the primordial plasma (T = 100 MeV) and bound on the 
neutrino mass (m Ur < 10 MeV). It is worth mention- 
ing that the decoupling temperature of r-neutrinos found 
in Ref. [2(| is p=> 3.7 MeV. Thus down to this tempera- 
ture r-neutrinos are in thermal equilibrium with primeval 
plasma. 

Although we studied the r-neutrino helicity relaxation 
in the early Universe, the interaction of the neutrino spin 
with the gravitational field (space-time curvature) was 
not taken into account in our work, i.e. all the calcu- 
lations were performed in the minkowskian space-time. 
However the neutrino propagation in strong gravitational 
fields of the specific configurations (e.g., Kerr metric) can 
cause the asymmetry between neutrinos and antineutri- 
nos. The most recent analysis of this issue is presented 
in Ref. [2]]]. In spite of the fact that the Universe is 
assumed to expand isotropically this effect is of great im- 
portance on the early stages of the Universe evolution 
when relic neutrinos interact with primordial black holes 
before neutrinos decoupled (see Ref. [2^1. This process 
could influence the neutrino helicity relaxation. 

It should be noted that for the first time the cosmo- 
logical upper bound on the sum of all neutrinos masses 
was derived in Ref. |23|. One can obtain the contribu- 
tion of stable light neutrinos to the matter density in the 
Universe and then compare the calculated value of the 
neutrino contribution with the observational value of to- 
tal matter density. The "classical" cosmological bound 
(Gerstein-Zeldovich limit) is very small 



< 14 eV. 



In our paper the production rate of " wrong-helicity" (i.e. 
right-handed) Dirac neutrinos was constrained by the 
Hubble parameter. For the first time the similar idea 
was used in Ref. [24| . It was assumed that if the popula- 
tion of the right-handed neutrino states in the primordial 
plasma is great, then they could influence the primordial 
nucleosynthesis by enlarging the effective number of the 
neutrino species. 

In conclusion we note that the helicity evolution of a 
neutrino interacting with matter having random distri- 
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butions of number density, velocities and polarizations 
has been studied in this paper. The basic equation was 
the generalized relativistic invariant Bargmann-Michel- 
Telegdi equation. We have derived the expression for the 
neutrino helicity relaxation rate which depends on the 
correlation functions of the hydrodynamical currents of 
background fermions. These correlators have been evalu- 
ated for the case of the ultrarelativistic electron-positron 
plasma. Then we have studied r-neutrino helicity evo- 
lution in plasma of the early Universe at T = 100 MeV. 
Comparing the helicity relaxation rate with the Hubble 
constant we have obtained the upper bound on the r- 



neutrino mass which turned out to be « 10 MeV. This 
r-neutrino mass limit is consistent with modern labora- 
tory constraints. 
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